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ABSTRACT 

The  Luzon  Strait  is  the  generation  region  for  strong  internal  tides  that  radiate  westward  into  the  South 
China  Sea  and  eastward  into  the  western  Pacific.  Intrusions  of  the  Kuroshio  and  strong  mesoscale  variability 
in  the  Luzon  Strait  can  influence  their  generation  and  propagation.  Here,  the  authors  use  eight  moorings  and 
two  numerical  models  to  investigate  these  relationships  by  quantifying  the  coherence  of  the  diurnal  and 
semidiurnal  internal  tides  in  the  Luzon  Strait.  This  study  finds  that  the  level  of  coherence  of  internal  tide 
generation,  energy,  and  energy  flux  is  quite  variable,  depending  on  the  specific  location  within  the  Luzon 
Strait.  Large-scale  spatial  patterns  in  internal  tide  pressure  and  velocity  exist  across  the  region,  shaped  by  the 
bathymetry,  mean  flow,  and  stratification.  Internal  tide  coherence  is  lower  (<30%)  near  large  gradients  in  this 
pattern  (predominantly  along  the  eastern  ridge),  which  are  shifted  by  the  variable  Kuroshio  and  mesoscale 
fields.  At  other  locations  within  the  Luzon  Strait,  the  internal  tide  is  largely  coherent  (>80%),  and  simple 
calculations  suggest  that  remote  sources  of  internal  tides  could  account  for  these  small  decreases  in  coherence. 
To  the  west  of  the  Luzon  Strait  (away  from  the  primary  generation  regions),  the  model  suggests  that  diurnal 
internal  tide  energy  is  more  coherent  than  semidiurnal. 


1.  Introduction 

Internal  tides,  generated  where  the  barotropic  tides 
force  stratified  water  over  underwater  bathymetry,  are 
thought  to  provide  a  significant  fraction  of  the  2  TW  of 
energy  required  to  maintain  abyssal  stratification  and 
the  meridional  overturning  circulation  (Munk  and 
Wunsch  1998).  They  contribute  to  ocean  mixing  through  a 
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cascade  of  energy  to  smaller-scale  waves  or  by  shoaling 
on  remote  continental  slopes.  A  better  understanding  of 
the  generation,  propagation,  and  dissipation  of  internal 
tides  near  major  generation  sites  is  necessary  to  un¬ 
derstand  and  predict  these  processes  under  present  and 
future  oceanographic  conditions. 

One  large  source  of  internal  tides  is  the  Luzon  Strait 
(Simmons  et  al.  2004;  Jan  et  al.  2008),  where  internal 
tides  are  generated  at  two  steep  submarine  ridges 
separating  the  South  China  Sea  (SCS)  from  the  Pacific 
Ocean  (Fig.  1).  In  addition  to  the  complex  bathymetry, 
the  Luzon  Strait  also  features  intrusions  of  the  Kur¬ 
oshio  and  high  mesoscale  variability.  The  Internal 
Waves  in  Straits  Experiment  (IWISE)  is  a  multi¬ 
investigator  initiative  funded  through  the  Office  of 
Naval  Research  to  investigate  internal  wave  dynamics 
in  this  complex  environment  using  a  combination  of 
shipboard  campaigns  and  moored  observations 
(Alford  et  al.  2015).  Two  major  scientific  objectives  of 
IWISE  are  determining  how  internal  tides  in  the  Lu¬ 
zon  Strait  are  affected  by  1)  a  double-ridge  system  with 
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Fig.  1.  (a)  Map  of  larger  region,  including  SCS.  White  box  indicates  region  shown  in  lower 
panels.  Black  diamonds  indicate  locations  of  the  eight  moorings  deployed  during  summer  2011. 
White  lines  through  moorings  indicate  the  northern  and  southern  lines  referred  to  in  the  text. 
Map  of  the  Luzon  Strait  with  mooring  locations  and  depth-integrated  (b)  diurnal  (Dl)  and 
(c)  semidiurnal  (D2)  energy  fluxes.  Modeled  MITgcm  energy  fluxes  are  shown  in  white.  Fluxes 
measured  from  LADCP/CTD  stations  are  shown  in  yellow.  Time-mean  fluxes  from  moorings 
A1  and  S9  are  shown  in  blue.  Boxes  around  these  moorings  indicate  regions  shown  in  Fig.  10. 
The  1000-m  isobath  is  contoured. 


complex  topography  and  2)  strong,  time-variable  cur¬ 
rents  and  stratification. 

As  it  propagates  west  from  the  ridges,  the  internal  tide 
often  steepens  into  nonlinear  internal  waves  (NLIWs) 
(Ramp  et  al.  2004;  Farmer  et  al.  2009;  Li  and  Farmer 
2011;  Lien  et  al.  2014)  that  are  associated  with  large 
vertical  displacements,  currents,  and  turbulence  that 
may  impact  the  local  ecosystem  (Wang  et  al.  2007; 
Moore  and  Lien  2007).  Improving  predictive  capability 
for  these  NLIWs  through  a  better  understanding  of  the 
variability  of  internal  tide  generation  at  the  Luzon  Strait 
is  another  motivation  for  this  work.  Park  and  Farmer 
(2013)  recently  showed  that  accounting  for  effects  of  the 
Kuroshio  on  wave  propagation  west  of  the  strait 


improves  the  prediction  of  NLIW  timing,  but  predicting 
the  amplitude  is  still  a  challenge. 

The  two  ridges  in  the  Luzon  Strait  contain  multiple 
internal  tide  generation  sites,  and  the  resultant  com¬ 
bination  of  waves  produces  interference  patterns  in 
depth-integrated  energy  fluxes  (Figs.  lb,c).  Energy 
flux  vectors  form  a  clockwise  gyre  in  the  northern 
section  of  the  strait,  where  both  the  west  and  east 
ridges  generate  significant  internal  tides.  To  the  south, 
where  the  western  ridge  is  deeper  and  does  not  gen¬ 
erate  strong  internal  tides,  energy  fluxes  emerge  pri¬ 
marily  from  the  east  ridge,  consistent  with  a  single 
source  there.  These  patterns  are  predicted  in  models 
(Jan  et  al.  2008;  Buijsman  et  al.  2012)  and  confirmed  in 
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measurements  made  during  the  IWISE  pilot  cruise 
(Alford  et  al.  2011). 

The  main  goal  of  this  work  is  to  quantify  influences 
from  the  mesoscale  or  Kuroshio  using  the  coherence  of 
internal  tides  over  several  spring-neap  cycles.  The  co¬ 
herent  portion  of  the  internal  tide  maintains  a  constant 
phase  and  amplitude  with  respect  to  the  barotropic  tidal 
forcing.  Quantifying  the  coherence  of  the  internal  tide  is 
important  for  several  reasons.  Estimates  of  internal  tide 
energy  flux  made  from  satellite  altimetry  only  capture 
the  coherent  portion  of  the  signal  (Ray  and  Zaron  2011; 
Zhao  2014).  Knowledge  of  the  coherent-incoherent 
partitioning  is  needed  to  make  accurate  estimates  of 
internal  tide  propagation  and  induced  mixing  or  to  re¬ 
move  the  internal  tide  signal  for  investigation  of  other 
nontidal  processes.  Coherence  is  also  closely  related  to  the 
predictability  of  internal  tides  and  associated  NLIWs.  If  it 
were  100%  coherent,  the  internal  tide  could  be  predicted 
given  knowledge  of  the  barotropic  forcing  and  the  (con¬ 
stant)  amplitude  and  phase  of  each  constituent.  Several 
mechanisms  could  lead  to  the  incoherence  of  the  internal 
tide  in  the  Luzon  Strait.  Although  they  are  not  necessarily 
independent,  we  will  here  divide  them  into  a  conceptual 
framework  of  “local”  and  “remote”  mechanisms. 

Local  refers  to  processes  occurring  in  the  Luzon  Strait, 
encompassing  both  ridges.  First,  generation  at  the  ridges 
may  be  modified  by  changes  in  stratification  or  shear  over 
those  sites.  The  Kuroshio  is  the  dominant  source  of  var¬ 
iability  in  stratification  and  shear  in  the  Luzon  Strait. 
Generation  at  one  ridge  may  also  be  modified  by  waves 
generated  at  the  other  ridge.  Changes  in  phase  speed 
between  the  ridges  caused  by  the  Kuroshio  could  alter  the 
phase  of  those  waves,  leading  to  variability  in  generation 
through  changes  in  the  bottom  perturbation  pressure 
(Kelly  and  Nash  2010;  Zilberman  et  al.  2011).  Finally,  it  is 
possible  that  the  3D  interference  pattern  in  the  northern 
portion  of  the  strait  (Buijsman  et  al.  2014)  could  be 
shifted  by  the  Kuroshio,  leading  to  variability  and  loss  of 
coherence  where  large  gradients  exist. 

Remote  refers  to  mechanisms  involving  signals  that 
come  from  outside  of  the  Luzon  Strait.  A  portion  of  the 
westward-propagating  diurnal  internal  tide  from  the 
Luzon  Strait  is  expected  to  reflect  back  from  the  conti¬ 
nental  shelf,  while  the  semidiurnal  internal  tide  is 
expected  to  propagate  onshelf  (Klymak  et  al.  2011). 
Mesoscale  variability  between  the  shelf  and  Luzon  Strait 
could  alter  the  phase  of  the  returning  signals,  leading  to 
lower  coherence.  Another  remote  source  is  the  Mariana 
Arc,  which  generates  strong  semidiurnal  internal  tides 
that  propagate  westward  to  the  Luzon  Strait.  These 
waves  propagate  through  a  region  of  strong  mesoscale 
variability  east  of  the  Luzon  Strait,  so  their  phase  upon 
arrival  could  be  highly  variable  (Kerry  et  al.  2013). 


In  this  paper,  we  characterize  the  internal  tide  struc¬ 
ture  and  variability  in  the  Luzon  Strait  over  several 
spring-neap  cycles,  focusing  on  two  moorings  with 
nearly  full-depth  coverage  over  50  days.  Data,  methods, 
and  oceanographic  context  are  presented  in  section  2. 
Results  are  presented  in  section  3,  focusing  on  detailed 
measurements  from  two  moorings.  In  section  4,  the  co¬ 
herent  and  incoherent  portions  of  the  internal  tide  are 
separated.  In  section  5,  we  examine  the  coherence  of  the 
internal  tide  in  the  model  and  estimate  the  effect  of 
several  mechanisms  on  the  internal  tide  coherence.  A 
summary  of  the  main  conclusions  is  given  in  section  6. 

2.  Data 

Observations  were  made  during  IWISE,  a  research 
initiative  to  study  the  generation,  propagation,  and  dis¬ 
sipation  of  internal  waves  in  the  Luzon  Strait.  A  pilot 
experiment  was  conducted  during  summer  2010,  con¬ 
sisting  of  two  short-term  moorings  and  19  shipboard 
LADCP/CTD  stations  (Alford  et  al.  2011)  along  a 
northern  and  southern  line.  During  the  summer  of  2011, 
an  eight-element  moored  array  (Fig.  1)  was  deployed 
along  the  same  northern  and  southern  lines.  The  goal  of 
the  array  was  to  occupy  the  same  locations  studied 
during  the  pilot  for  a  longer  period  in  order  to  measure 
the  time  variability  of  those  signals  over  several  spring- 
neap  cycles.  As  a  result  of  the  challenging  environment, 
not  all  moorings  returned  full-depth  data.  In  this  paper, 
we  focus  on  two  moorings  (Al  and  S9)  where  nearly  full- 
depth  measurements  were  obtained,  allowing  reliable 
separation  of  barotropic  and  baroclinic  signals  and  cal¬ 
culation  of  internal  tide  energy  and  flux. 

a.  Al 

Mooring  Al  was  located  on  the  eastern  ridge  in  the 
central  Luzon  Strait  (20°31.4'N,  121o40.5'E)  at  a  water 
depth  of  1446  m.  A  “Stablemoor”  platform  at  a  depth  of 
705  m  housed  upward-  and  downward-looking,  75-kHz 
RDI  Long  Ranger  acoustic  Doppler  current  profilers 
(AD CPs),  which  measured  velocity  over  depths  of 
—120-650  and  —725-1200  m  with  5-min  time  resolution. 
An  upward-looking  300-kHz  ADCP  at  80  m  measured 
near-surface  velocity.  Seabird  (SBE)-37  MicroCATs 
(temperature,  salinity,  and  pressure)  and  SBE-56s 
(temperature)  were  deployed  along  the  mooring  with 
2-min  sampling  rate.  Pressure  records  from  the  SBE-37 
MicroCATs  were  used  to  correct  for  mooring 
knockdown. 

b.  S9 

Mooring  S9  was  located  on  the  southern  portion  of  the 
eastern  ridge  (19°18'N,  121°1.8/E)  at  a  water  depth  of 
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2286  m.  Two  stacked  McLane  moored  profilers  (MPs) 
profiled  depth  ranges  of  330-1260  and  1300-2220  m, 
measuring  horizontal  velocity,  temperature,  and  salin¬ 
ity.  The  MP  profiles  at  a  nominal  speed  of  0.25 ms-1, 
resulting  in  a  complete  profile  approximately  every 
hour.  The  MP  has  a  vertical  resolution  of  —10  m  for 
velocity  and  2  m  for  density.  A  subsurface  float  at  300  m 
contained  an  upward-looking,  75-kHz  ADCP  that 
measured  velocity  from  —30  to  300  m.  Between  this  float 
and  a  smaller  float  at  70  m,  a  T  chain  consisting  of  SBE- 
56  temperature  loggers  spaced  15  m  apart  measured 
temperature.  Pressure  records  from  the  bottom  and  top 
of  the  T  chain  were  used  to  correct  the  SBE-56  tem¬ 
perature  records  for  mooring  knockdown  (see  the  ap¬ 
pendix  for  details).  Salinity  was  estimated  for  the  T 
chain  using  the  observed  T-S  relation  from  the  SBE-37. 
For  the  following  analysis,  all  data  from  both  moorings 
were  interpolated  onto  a  uniform  grid  with  10-m  vertical 
and  1.5-h  temporal  spacing. 

c.  Models 

We  also  make  use  of  two  numerical  models  in  this 
analysis.  The  first  is  the  Massachusetts  Institute  of 
Technology  general  circulation  model  (MITgcm) 
(Marshall  et  al.  1997),  which  uses  more  detailed 
bathymetry  and  provides  higher  spatial  resolution  but 
does  not  contain  any  mesoscale  variability.  The  second 
is  the  Luzon  Strait  Ocean  Nowcast/Forecast  System 
(LZSNFS,  hereinafter  LZS)  model,  which  has  coarser 
resolution  but  is  run  over  a  longer  time  period  and 
contains  a  realistic  Kuroshio  and  mesoscale  variability 
in  addition  to  tidal  forcing. 

1)  MITgcm  model 

MITgcm  was  used  in  3D  configuration,  as  described  in 
Buijsman  et  al.  (2014).  The  model  uses  realistic  topog¬ 
raphy  merged  from  high-resolution  gridded  multibeam 
data  with  a  resolution  of  — 300-m  and  30-arc-s  resolution 
global  topography/bathymetry  grid  (SRTM30_PLUS) 
data  from  the  Smith  and  Sandwell  database  with  a  res¬ 
olution  of  —1  km  (Smith  and  Sandwell  1997).  The  den¬ 
sity  stratification  is  horizontally  uniform  and  derived 
from  temperature  and  salinity  data  collected  in  between 
the  ridges  (Alford  et  al.  2011).  The  tidal  forcing  at 
the  east,  west,  north,  and  south  model  boundaries  is 
composed  of  barotropic  velocities  extracted  from 
TOPEX/Poseidon,  version  7.2  (TPX07.2;  Egbert  and 
Erofeeva  2002). 

2)  LZS  MODEL 

The  LZS  model  is  an  application  of  the  NRL  Ocean 
Nowcast/Forecast  System  (Ko  et  al.  2008).  The  system 
is  an  integration  of  a  dynamical  ocean  model  and  a 


statistical  data  analysis  model.  The  model  derives  its 
ocean  bottom  topography  from  the  Digital  Bathyme¬ 
try  Data  Base  with  2-min  resolution.  Open  boundary 
conditions  are  derived  from  a  larger-scale  model  for 
the  entirety  of  the  East  Asian  seas  (Ko  et  al.  2009). 
The  barotropic  tidal  currents  of  eight  tidal  constitu¬ 
ents  (Ki,  Oi,  Pi,  Qi,  K2,  M2,  N2,  and  S2)  are  obtained 
from  TPX07.2  (Egbert  and  Erofeeva  2002).  This 
model  includes  the  reflection  of  the  internal  tide  from 
the  continental  shelf  to  the  west  but  does  not  contain 
remote  internal  tides  from  outside  the  domain,  in¬ 
cluding  the  Mariana  Arc.  The  LZS  model  has  now 
been  used  in  several  studies  of  internal  waves  in  the 
Luzon  Strait  and  South  China  Sea  (Chen  et  al.  2013; 
Ma  et  al.  2013). 


3.  Results 

Our  primary  focus  is  the  internal  tide  energy  and  en¬ 
ergy  flux.  However,  we  first  take  a  moment  to  present 
the  basic  data  (Fig.  2)  and  velocity  spectra  (Fig.  3). 

a.  Basic  signals 

The  nearly  full-depth  coverage  allows  us  to  un¬ 
ambiguously  separate  the  barotropic  and  baroclinic 
signals.  Barotropic  tidal  velocity  is  defined  as 


UBT 


1 

H 


[u  (z,  t )  dz 

-H 


u(z)]  dz. 


(1) 


where  u(z,  t)  is  the  instantaneous  velocity,  and  u (z)  is  a 
slowly  varying  sub  tidal  component.  Baroclinic  velocity 
is  defined  as 


u'  =  u(z,f)  -UBT.  (2) 

The  spring-neap  cycle  of  tidal  forcing  is  evident  in  the 
barotropic  velocity  (Fig.  2a).  The  barotropic  velocities 
are  very  strong,  up  to  0.48  m  s-1  at  Al,  and  0.31  m  s-1  at 
S9.  Total  velocities  (Figs.  2b-d)  are  dominated  by  tidal 
frequency  motions  and  a  strong  spring-neap  cycle  as 
well  as  subtidal  motions.  Vertical  displacements  are 
tidally  dominated  and  increase  with  depth,  as  indicated 
by  the  variability  of  density  contours. 

b.  Spectra 

Rotary  spectra  of  baroclinic  horizontal  velocity  at  a 
depth  of  250  m  from  S9  are  presented  in  Fig.  3  (spectra 
from  Al  are  similar  in  character  and  are  not  shown).  The 
spectra  are  dominated  by  peaks  at  tidal  frequencies  as 
well  as  sum  and  difference  frequencies  (e.g.,  M2  +  Ki) 
and  higher  harmonics.  The  spectra  are  elevated  about 
10-15  times  above  the  theoretical  Garrett-Munk 
(GM76)  internal  wave  spectrum  (Garrett  and  Munk 
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Fig.  2.  Raw  data  from  moorings  A1  and  S9.  (a)  Zonal  component  of  barotropic  velocity  [Eq.  (1)]. 
(b)  A1  zonal  velocity,  (c)  A1  meridional  velocity,  (d)  S9  zonal  velocity,  and  (e)  S9  meridional 
velocity.  Black  lines  are  isopycnal  contours  at  approximately  500-m  intervals.  Contour  values 
(kgnT3)  are  (1023.2, 1026.8,  and  1027.4)  for  A1  and  (1023.1, 1027.0, 1027.45,  and  1027.6)  for  S9. 


1972;  Cairns  and  Williams  1976)  but  follow  a  similar 
slope  and  fall  off  with  increasing  frequency.  The  extreme 
elevation  above  GM  levels  nearly  completely  masks  the 
usually  prominent  inertial  peak,  which  here  is  only 
barely  detectable. 

c.  Mesoscale  velocity  and  Kuroshio 

We  use  the  large  spatial  extent  of  the  mooring  array 
and  the  LZS  model  to  characterize  the  spatial  and 
temporal  structure  of  the  Kuroshio  and  mesoscale  ve¬ 
locity  field  during  the  experiment.  Figure  4  shows  low- 
passed  (cutoff  of  0.25  cpd)  velocity  vectors  from  moored 
ADCPs  (color)  at  150  m,  with  the  color  indicating  the 
time  of  each  measurement  (only  every  3  days  are  shown 
for  clarity).  Gray  vectors  are  the  time-mean  low-passed 
velocity  at  150  m  from  the  LZS  model  during  the 


experiment,  which  are  generally  consistent  with  the 
low-pass  filtered  observations.  A  persistent  Kuroshio 
intrusion  was  present  during  the  experiment,  entering 
the  Luzon  Strait  from  the  southeast  and  passing  over 
the  western  half  of  the  strait  near  the  northern  line 
before  exiting  to  the  east  of  Taiwan.  It  can  be  seen  in 
observations  at  moorings  N2  and  MPN,  where  low- 
frequency  velocity  was  strong  and  to  the  north  during 
the  experiment.  Note  that  A1  was  not  in  the  path  of  the 
Kuroshio  intrusion,  which  appears  to  be  deflected 
around  the  island  and  shallow  topography  to  the 
southeast  of  the  mooring.  Velocity  measured  at 
moorings  along  the  southern  line  was  also  outside  of 
the  main  Kuroshio  intrusion  but  exhibited  strong, 
variable  low-frequency  flows,  likely  eddies  spun  off  of 
the  Kuroshio  intrusion. 
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Fig.  3.  Rotary  spectrum  of  baroclinic  velocity  at  250  m  at 
mooring  S9.  Spectra  of  the  raw  velocity  measured  by  an  upward- 
looking,  75-kHz  ADCP  are  shown  in  black  (clockwise)  and  gray 
(anticlockwise).  The  clockwise  spectrum  of  the  gridded  MP  ve¬ 
locity  is  plotted  in  blue.  Gray  shading  shows  the  (scaled)  frequency 
response  of  the  bandpass  filter  used  to  isolate  diurnal  and  semi¬ 
diurnal  frequencies.  Red  line  is  15  X  the  theoretical  GM76  spec¬ 
trum  for  the  mean  N  at  that  depth.  Black  lines  at  bottom  of  figure 
indicate  the  95%  confidence  limits. 
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Fig.  4.  Low-frequency  velocity  across  the  Luzon  Strait  during 
the  experiment.  Colored  vectors  show  low-pass  filtered  velocities 
at  150  m  from  the  mooring  array.  Gray  vectors  are  time-mean 
velocity  vectors  at  150  m  from  the  LZS  model  during  the 
experiment. 


Low-frequency  velocity  and  stratification  at  S9  and 
A1  are  computed  by  low-pass  filtering  with  a  cutoff  of 
0.25  cpd  (Fig.  5).  Both  moorings  have  a  stratification 
maximum  near  150  m.  In  general  N 2,  defined  as 
(- g/p0)(dp/dz ),  where  g  is  the  acceleration  due  to 
gravity,  p  is  potential  density,  and  pQ  is  a  reference 
density,  is  higher  at  S9  than  at  Al.  Variable  N 2  is  seen  in 
the  upper  200-300  m  at  both  moorings.  The  low- 
frequency  velocities  are  surface  intensified  at  both 
moorings.  The  velocity  at  Al  appears  to  have  a  weak 
spring-neap  cycle,  likely  because  of  the  tidal  rectifica¬ 
tion.  The  low-frequency  velocity  at  S9  is  dominated  by 
strong  (0.5 ms-1),  near-surface  (0-300 m)  pulses  with 
time  scales  of  a  few  days. 

d.  Internal  tides 

Internal  tide  energy  and  energy  flux  are  computed 
using  standard  methods  (Kunze  et  al.  2002;  Nash  et  al. 
2005).  The  diurnal  and  semidiurnal  frequency  bands  are 
separated  by  applying  a  bandpass  filter  (gray  in  Fig.  3)  to 
the  velocity  and  displacement  data.  The  shallowest  and 
deepest  velocities  are  extended  by  linear  extrapolation 
to  the  surface  and  bottom,  respectively.  Gaps  in  velocity 
were  limited  to  the  upper  20  m  and  lower  100  m  at  Al, 
and  the  upper  40  m  and  lower  50  m  at  S9.  Stratification 
N 2  and  vertical  displacement  17  are  linearly  interpolated 
from  the  shallowest  (deepest)  measured  value  to  zero  at 
the  surface  (bottom).  Gaps  in  displacement  were  limited 


to  the  upper  70  m  and  lower  45  m  at  Al  and  the  upper 
90  m  and  lower  60  m  at  S9.  Smaller  gaps  in  velocity  and 
displacement  in  the  middle  of  the  water  column  are  filled 
by  linear  interpolation. 

Horizontal  kinetic  energy  (HKE)  and  available  po¬ 
tential  energy  (APE)  are  computed  as 


HKE  =  f  p(m/2  4-  v12),  and 

(3) 

APE  =  ^pV2(tj2)  , 

(4) 

where  17  is  vertical  displacement  computed  over  3-day 
windows  and  brackets  indicate  averaging  over  a 
wave  period. 

The  perturbation  pressure  p'  is  computed  as 

f° 

P'=^surf+  J  P'gdZ, 

(5) 

where  the  density  perturbation  is  given  as 

P  =  -  N2 7]. 

8 

(6) 

The  surface  pressure  psurf  is  not  measured  but  can  be 
inferred  from  the  baroclinicity  condition  that  the  depth- 
averaged  pressure  perturbation  must  vanish  (Althaus 
et  al.  2003): 
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Fig.  5.  Low-pass  filtered  stratification  and  velocity  (4-day  cutoff)  at  S9  and  Al.  (a)  Zonal 
barotropic  velocity  at  S9  (red)  and  Al  (black),  (b)  Stratification  at  Al.  (c)  Stratification  at  S9. 
(d)  Zonal  velocity  at  Al.  (e)  Meridional  velocity  at  Al.  (f)  Zonal  velocity  at  S9.  (g)  Meridional 
velocity  at  S9. 


J  HP'(Z’t)dz  =  °-  (7) 

Determining  psurf  and  p'  requires  integrating  over  the 
full  water  depth.  Thus,  vertical  gaps  in  measurements 
can  introduce  errors  in  this  quantity  (Nash  et  al.  2005; 
Rainville  and  Pinkel  2006).  Sensitivity  studies  showed 
that  errors  in  psurf  and  p'  because  of  the  interpolating 
surface  gaps  were  small  (<10%),  allowing  us  to  compute 
accurate  profiles  of  p'  over  the  depth  ranges  we  sampled. 

Energy  flux  is  defined  as  the  covariance  of  velocity 
and  pressure  perturbations: 

F=(uV).  (8) 

Although  profiles  of  p'  over  the  measured  depth  range 
are  accurate,  the  depth-integrated  energy  and  flux 


magnitude  may  differ  from  the  true  value  because  of  the 
vertical  gaps.  This  difference  depends  on  the  depth 
profile  of  the  energy  flux,  which  is  normally  surface  in¬ 
tensified.  To  estimate  this  error,  we  resampled  MITgcm 
model  data  from  the  mooring  locations  and  computed 
energy  flux  in  the  same  manner  as  the  observations. 
Errors  from  the  true  value  were  1%  and  17%  for  diurnal 
and  semidiurnal  at  Al  and  5%  and  9%  for  diurnal  and 
semidiurnal  fluxes  at  S9.  Importantly,  these  errors  give 
rise  to  largely  time-independent  biases  that  do  not  affect 
our  present  focus  on  the  time-variable  signals  and  their 
coherence. 

We  also  compute  barotropic  to  baroclinic  conversion, 
defined  as  (Kelly  and  Nash  2010) 

C={\H-(uBTp'b)),  (9) 
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Fig.  6.  Diurnal  internal  tide  quantities  at  S9:  (a)  barotropic  velocity  uBT,  (b)  baroclinic  zonal  velocity, 
(c)  perturbation  pressure,  (d)  energy  (HKE  +  APE),  (e)  zonal  energy  flux,  and  (f)  meridional  energy  flux. 


where  \H  is  the  bottom  slope,  p'b  is  the  perturbation 
pressure  at  the  bottom,  and  brackets  indicate  averaging 
over  a  wave  period. 

e.  Depth-time  structure  of  internal  tides 

We  first  examine  and  contrast  the  depth-time 
structure  of  the  diurnal  internal  tide  at  S9  and  Al. 
Semidiurnal  quantities  are  omitted  for  brevity;  see  the 
online  supplemental  material  for  figures.  The  diurnal 
signals  are  noticeably  more  coherent  at  S9,  which  we 
will  further  quantify  in  a  future  section.  The  diurnal 


internal  tide  at  S9  (Fig.  6)  exhibits  a  relatively  clean 
spring-neap  cycle  because  of  the  beating  of  the  K4  and 
Oi  components.  The  vertical  structure  of  perturbation 
pressure,  energy,  and  flux  appears  to  be  low  mode  and 
is  similar  during  each  spring  tide.  Energy  and  flux  are 
largest  in  the  upper  500  m,  with  a  smaller  peak  in  the 
lower  water  column. 

The  diurnal  velocity  at  Al  (Fig.  7b)  also  shows  a 
strong  spring-neap  cycle.  Energy  (dominated  here  by 
HKE)  has  a  similar  vertical  structure  during  all  four 
spring  tides.  However,  perturbation  pressure  (c)  is  more 
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Fig.  7.  As  in  Fig.  6,  but  for  Al. 
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variable  in  both  time  and  depth  and  does  not  have  such 
an  obvious  spring-neap  cycle.  As  a  result,  the  depth  and 
time  structure  of  flux  is  more  variable,  notably  during 
the  last  spring  tide  (yearday  >  205).  It  will  be  shown  in  a 
later  section  that  this  lowers  the  coherence  of  the  diurnal 
internal  tide  at  Al. 

f  Modal  decomposition 

The  modal  structure  of  the  internal  tide  was  also  ex¬ 
amined  by  fitting  velocity  and  displacement  profiles  to 
mode  shapes  determined  numerically  from  the  IW  ver¬ 
tical  structure  equation: 


^<Kz)  +  ^4(z)  =  0,  (10) 

dzz  cln 

where  4>  is  the  vertical  mode  shape,  and  cn  is  the 
eigenspeed.  The  modal  distributions  of  time-averaged, 
depth-integrated  APE,  HKE,  and  energy  flux  at  S9  and 
Al  are  shown  in  Fig.  8.  The  modal  distributions  are  red, 
with  most  of  the  variance  contained  in  the  first  three 
modes.  Energy  flux  is  dominated  by  mode  1  with  the 
exception  of  the  diurnal  internal  tide  at  Al,  where  mode 
1  is  the  largest  but  modes  2  and  3  make  up  a  larger 
fraction  of  the  total  (as  can  be  seen  by  the  cumulative 
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Fig.  8.  Modal  distribution  of  (top)  APE,  (middle)  HKE,  and  (bottom)  energy  flux  at  (left)  S9 
and  (right)  Al.  Lines  with  circles  show  the  cumulative  percent  total,  indicated  by  horizontal 
bars  on  right  of  axes  and  scale  at  upper  right. 


percent  totals).  The  dominance  of  low-mode  flux  is  typical 
of  internal  tides,  since  the  pressure  tends  to  act  as  a  low- 
pass  filter  (Nash  et  al.  2005).  HKE  at  both  moorings  has 
more  variance  at  higher  modes  than  energy  flux  but  still 
tends  to  have  a  red  spectrum,  with  about  75%  of  the  total 
HKE  contained  in  the  first  2  modes.  At  S9,  the  mode-1 
APE  contains  about  75%  of  the  total  in  the  diurnal  band 
and  less  than  50%  in  the  semidiurnal  band.  At  Al, 
semidiurnal  APE  is  dominated  by  mode-1  but  diurnal 
APE  is  distributed  more  evenly  among  modes  (note 
though  that  diurnal  energy  at  Al  is  dominated  by  HKE). 

g.  Depth-integrated  internal  tide  energy  and  flux 

Next,  we  examine  time  series  of  depth-integrated  in¬ 
ternal  tide  energy,  energy  flux,  and  conversion  (Fig.  9). 
Time-mean  quantities  are  summarized  in  Table  1.  At  S9, 
the  diurnal  barotropic  forcing  exceeds  the  semidiurnal 
forcing  by  about  a  factor  of  2.  This  is  reflected  in  the 
greater  diurnal  magnitude  of  energy  (Fig.  9b),  flux 
(Fig.  9c),  and  conversion  (Fig.  9d)  at  S9. 

At  Al  (gray),  diurnal  and  semidiurnal  barotropic 
forcing  have  comparable  magnitudes.  Depth-integrated 


energy  is  large  in  both  bands,  with  semidiurnal  exceed¬ 
ing  diurnal  most  of  the  time.  Note  that  the  diurnal  en¬ 
ergy  at  Al  is  dominated  by  HKE.  While  the  diurnal  and 
semidiurnal  barotropic  forcing  and  energy  magnitudes 
are  comparable,  the  energy  flux  and  conversion  are 
dominated  by  the  semidiurnal  band.  The  diurnal  energy 
flux  is  small  compared  to  the  semidiurnal  flux  but  still 
very  large  compared  to  typical  open-ocean  values. 
Semidiurnal  flux  appears  to  be  well  correlated  with  the 
barotropic  forcing.  Diurnal  flux  appears  less  correlated 
with  barotropic  forcing,  even  though  the  diurnal  energy 
is  large  and  well  correlated. 

A  map  view  of  the  depth-integrated  diurnal  (left)  and 
semidiurnal  (right)  energy  fluxes  at  Al  (top)  and  S9 
(bottom)  is  shown  in  Fig.  10.  As  also  seen  in  Fig.  1,  the 
observed  and  modeled  fluxes  generally  show  the  same 
patterns,  with  an  interference  pattern  to  the  north 
(Figs.  10a, b)  and  a  single  source  to  the  south  (Figs.  10c, d). 
The  observed  diurnal  fluxes  at  Al  (Fig.  10a)  show  the 
greatest  directional  variability,  with  some  of  the  fluxes 
toward  the  southwest  in  agreement  with  the  model  but 
also  directed  west  and  east  at  times.  The  semidiurnal 
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Fig.  9.  Time  series  of  depth-integrated  internal  tide  quantities  at  moorings  S9  (black)  and  A1 
(gray).  Diurnal  quantities  are  shown  with  solid  lines  and  semidiurnal  are  shown  with  dashed 
lines,  (a)  Baro tropic  velocity  amplitude,  (b)  Depth-integrated  energy  (HKE  +  APE),  (c)  Depth- 
integrated  energy  flux  magnitude,  (d)  Barotropic  to  baroclinic  conversion. 


band  at  S9  also  shows  more  directional  variability  but 
less  than  the  diurnal  band  at  Al. 

4.  Coherent/incoherent  separation 

We  next  decompose  the  internal  tide  signals  into  co¬ 
herent  and  incoherent  components  using  a  method 


similar  to  Nash  et  al.  (2012).  A  harmonic  fit  to  the  major 
diurnal  (Ki,  Oi,  and  Qi)  or  semidiurnal  (M2,  S2,  and  N2) 
tidal  constituents  is  performed  on  the  bandpass-filtered 
time  series  of  diurnal  or  semidiurnal  velocity  and  pres¬ 
sure  at  each  depth.  Note  that  the  bandwidth  of  the  filter 
will  affect  the  partitioning  of  coherent  and  incoherent 


Table  1.  Time-mean,  depth-integrated  energy,  and  energy  flux  for  diurnal  and  semidiurnal  internal  tide  at  moorings  S9  and  Al.  Direction 

is  in  degrees  true  (clockwise  from  north). 


Diurnal 

Semidiurnal 

E  (kJ  m-2) 

]ff  (kW  m-1) 

d  (°True) 

E  (kJ  m-2) 

jf[  (kW  m-1) 

6  (True) 

S9 

24 

25 

339 

15 

14 

313 

Al 

21 

5 

236 

27 

27 

169 
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Fig.  10.  (left)  Diurnal  and  (right)  semidiurnal  depth-integrated  energy  fluxes  at  moorings 
(top)  A1  and  (bottom)  S9.  White  vectors  are  time-mean  modeled  (MITgcm)  fluxes.  Observed 
fluxes  are  colored  by  time. 


signals  (since  an  infinitely  thin  passband  contains  just 
one  frequency  and  is  100%  coherent  by  definition).  The 
passbands  we  used  (shown  in  gray  in  Fig.  3)  are  0.85-1.15 
cpd  for  the  diurnal  band  and  1. 7-2.3  cpd  for  the  semi¬ 
diurnal  band,  which  include  the  three  diurnal  or  semi¬ 
diurnal  constituents  used  for  the  harmonic  fits.  We  do 
not  include  Pi  or  K2  in  the  fits  because  the  record  length 
is  not  long  enough  to  resolve  them  from  or  S2,  re¬ 
spectively.  The  coherent  part  ( u ,  p)c  is  defined  as  the 
reconstructed  fit,  and  the  remainder  is  defined  as  the 
incoherent  part: 


Velocities  are  separated  into  barotropic  and  baroclinic 
components: 

u'c  =  uc  —  and  (14) 

u'i  =  ut  —  ufT ,  (15) 

where  uf  j  is  the  depth  mean  of  uCj.  Then  baroclinic 
energy,  energy  flux,  and  conversion  are  computed  sep¬ 
arately  for  each  component.  In  addition  to  the  coherent 
and  incoherent  terms,  there  is  a  third  component  re¬ 
sulting  from  the  cross  terms.  APE  is  given  by 


m  .  =  u  —  uc, 
v  •  =  v  —  vc,  and 

Pi=P-Pc • 


(11) 

(12) 

(13) 


APE  =  APE,  +  APE.  +  APE, 


APE  =  —  -  o' 2 
<  2  PN2  Pc  ’ 


(16) 

(17) 
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Fig.  11.  Results  of  coherent/incoherent  separation  for  diurnal  band  at  S9.  (a)  Zonal 
barotropic  velocity,  (b)  HKE.  (c)  APE.  (d)  Zonal  energy  flux,  (e)  Meridional  energy  flux, 
(f)  Conversion.  Coherent,  incoherent,  and  cross-term  components  are  shown  in  color. 
Dashed  black  line  shows  the  total  of  each  quantity.  Number  in  upper  left  of  each  panel  indicates 
the  percent  of  variance  explained  by  coherent  portion  [Eq.  (32)].  Note  that  the  y-axis  limits  are 
not  equal  in  Figs.  11-14. 


“d 

p2 

(18) 

HKEcr  =  piu^Uj  +  v'cv'j) . 

Energy  flux  is  given  by 

(23) 

APE 

(19) 

L  =  C  +  ^  +  Fcr, 

Fc  =  (ucPc)  > 

(24) 

(25) 

HKE  is  given  by 

Fi  =  ( u'tp'i ),  and 

(26) 

HKE  =  HKEc  +  HKE;  +  HKEcr , 

(20) 

51 

s 

+ 

II 

(27) 

HKEC  =  1  p(uf  +  vf), 

(21) 

with  similar  expressions  for  the  meridional  component 
Fv.  Barotropic  to  baroclinic  conversion  is  given  by 

HKE(.  =  1  p{uf+vf),  and 

(22) 

C-Cc  +  Ci  +  C cr> 

(28) 
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Fig.  12.  As  in  Fig.  11,  but  for  the  semidiurnal  band  at  S9. 


II 

■fa 

H 

< 

$ 

(29) 

c.  =  ufT  ■  VHp'i,  and 

(30) 

ccr  =  UfT  ■  VHp'i  +  UfT  ■  \Hp'c . 

(31) 

Results 

As  expected,  because  of  their  proximity  to  strong 
internal  tide  generation  sites,  the  coherent  term 
(Figs.  11-14,  orange)  is  generally  the  largest  of  the 
three  components.  Although  not  strictly  tied  to  the 
barotropic  forcing,  we  loosely  interpret  the  coherent 
component  as  the  locally  generated  internal  tide,  in 
the  absence  of  changes  in  mesoscale  stratification  or 
velocity.  The  magnitude  and  direction  of  the  coherent 
flux  vary  according  to  the  spring-neap  cycle,  and 
spring  tides  are  not  identical  because  of  the  Qi  and  N2 
components. 


The  remaining  incoherent  and  cross  terms  are  deviations 
from  this  regular  cycle.  They  tend  to  also  be  the  largest 
during  spring  tides,  but  their  phase  and/or  amplitude  vary 
on  time  scales  shorter  than  the  ~  50-day  record,  and  as  a 
result  they  are  not  captured  in  the  harmonic  fits.  Fits  over 
individual  spring  tides  (not  shown)  are  more  coherent, 
suggesting  that  changes  in  the  internal  tide’s  amplitude  or 
phase  occur  on  time  scales  between  14  and  50  days. 

To  compare  the  coherence  at  different  locations,  we 
compute  the  percent  variance  of  each  quantity  that  is 
explained  by  the  coherent  portion  of  the  signal: 


%varx 


var(X-XJ 

var(X) 


X100, 


(32) 


which  is  indicated  in  the  upper  left  of  each  panel  in 
Figs.  11-14  and  Table  2. 

The  coherence  of  the  diurnal  (Dl)  energy  flux  at  A1  is 
less  than  30%,  significantly  less  than  the  diurnal  internal 
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Fig.  13.  As  in  Fig.  11,  but  for  the  diurnal  band  at  Al. 


tide  at  S9  and  the  semidiurnal  (D2)  internal  tide  at  both 
locations.  The  diurnal  Al  energy  flux  (Figs.  13d, e)  is 
dominated  by  the  cross-term  fluxes,  which  are  of  compa¬ 
rable  magnitude  as  the  coherent  fluxes  and  change  the  sign 
of  both  components  of  energy  flux,  resulting  in  highly 
variable  energy  flux  vectors  (as  seen  in  Fig.  10).  Further 
decomposing  the  cross-term  fluxes  into  their  two  compo¬ 
nents  [Eq.  (27)]  reveals  that  they  are  dominated  by  the  term 
involving  incoherent  pressure  (u'jj'j) .  This  is  consistent  with 
the  qualitative  incoherent  nature  of  perturbation  pressure 
and  energy  flux  seen  in  Fig.  7c.  Conversion  (Fig.  13f),  which 
depends  on  the  magnitude  and  phase  of  bottom  per¬ 
turbation  pressure,  is  also  less  coherent  at  Al  (63%). 

Cross-term  fluxes  in  the  semidiurnal  band  at  Al  and 
both  bands  at  S9  are  of  similar  magnitude  but  smaller 
relative  to  coherent  fluxes.  At  S9,  the  diurnal  zonal  Fcr 
(Fig.  lid)  is  negative  during  the  first  two  springs  and 
positive  during  the  last  two  springs.  Meridional  Fcr  is 
mostly  positive  during  the  first  two  springs,  and  negative 


during  the  last  two  springs.  The  zonal  Fcr  appears  to  be 
somewhat  correlated  with  the  mean  low-frequency  ve¬ 
locity  over  the  upper  1000  m  (Fig.  5e),  suggesting  that 
the  internal  tide  energy  flux  may  be  modified  by  the  low- 
frequency  eddy/mesoscale  velocities. 

The  levels  of  coherence  are  summarized  in  Fig.  15,  in 
which  it  can  be  seen  that  the  diurnal  band  at  S9  and  the 
semidiurnal  band  at  Al  are  largely  coherent  (>70%  for 
all  quantities).  Semidiurnal  motions  at  S9  are  slightly 
less  coherent  (>60%).  Diurnal  motions  at  Al  stand  out 
as  significantly  less  coherent  (<30%,  except  for  HKE) 
than  the  other  three,  consistent  with  the  visual  appear¬ 
ance  of  pressure  (Fig.  7c)  and  the  obvious  veering  of  the 
flux  direction  (Fig.  10a).  We  discuss  reasons  for  this  next. 

5.  Discussion 

Which  of  the  several  possible  physical  mechanisms  are 
responsible  for  the  observed  levels  of  internal  tide 
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Fig.  14.  As  in  Fig.  13,  but  for  the  semidiurnal  band  at  Al. 


incoherence  in  the  Luzon  Strait?  We  first  use  a  simple 
2D  model  to  estimate  the  effect  variable  stratification 
and  velocity  due  to  the  Kuroshio  has  on  internal  tides 
propagating  between  the  ridges  at  20.5°N.  The  model 
consists  of  a  coherent  mode-1  internal  tide  generated  at 
both  ridges,  with  the  total  observed  internal  tide  at  Al 
then  consisting  of  the  superposition  of  the  locally  gen¬ 
erated  wave  (the  observed  coherent  internal  tide)  and 
the  wave  arriving  from  the  western  ridge  after  traveling 
through  variable  phase  speed  between  the  ridges.  The 
phase  speed  is  calculated  from  the  Taylor-Goldstein 
equation  using  low-passed  velocity  and  stratification 
from  the  LZS  model.  The  resulting  coherence  of  flux 
and  conversion  in  both  bands  [Eq.  (32)]  is  99%,  showing 
that  this  mechanism  has  a  negligible  effect. 

Next,  we  use  a  similar  model  to  consider  the  impact 
of  the  diurnal  internal  tide  reflected  from  the  conti¬ 
nental  shelf  to  the  west.  In  this  scenario,  a  local 


(coherent)  internal  tide  is  generated  at  Al  and  prop¬ 
agates  west  to  the  slope  at  117.5°E,  where  V3  of  the 
energy  reflects  back,  consistent  with  Klymak  et  al. 
(2011).  The  resulting  diurnal  coherence  in  this  2D 
model  is  86%  for  flux  and  99%  for  conversion.  This 
mechanism  could  explain  some  of  the  lower  diurnal 
coherence  in  flux  at  Al,  but  the  observed  coherence  is 
significantly  lower.  However,  it  is  enough  to  explain 
the  88%  coherence  of  the  diurnal  flux  at  S9.  Klymak 


Table  2.  Percent  variance  of  observed  mode-1  energy,  flux,  and 
conversion  explained  by  coherent  portion  at  moorings  Al  and  S9. 


APE 

HKE 

Fu 

Fv 

F 

C 

Al-Dl 

16 

80 

38 

30 

24 

31 

A1-D2 

83 

79 

53 

85 

84 

88 

S9-D1 

90 

72 

85 

84 

88 

94 

S9-D2 

65 

60 

57 

61 

74 

67 
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Fig.  15.  Percent  coherence  of  internal  tides  at  A1  and  S9  in  both 
frequency  bands,  computed  via  Eq.  (32). 


et  al.  (2011)  note  that  the  amount  of  diurnal  energy 
reflected  is  sensitive  to  the  phase  between  mode-1 
and  mode-2  waves.  Using  the  LZS  model  data,  we 
calculated  the  mode-1  and  mode-2  phase  speeds,  and 
the  phase  difference  between  them  after  propagat¬ 
ing  to  the  shelf  (assuming  they  start  in  phase).  The 
phase  difference  has  a  range  of  about  80°,  suggesting 
that  the  fraction  of  energy  reflected  from  the  shelf 
could  be  quite  variable  in  time,  helping  to  explain 
the  modest  incoherence  observed  in  the  S9  diurnal 
internal  tide. 

Finally,  we  consider  the  effect  of  remote  (semi¬ 
diurnal)  waves  from  the  Mariana  Arc.  In  recent  mod¬ 
eling  studies,  Kerry  et  al.  (2014)  found  that  including 
remote  waves  from  the  Mariana  Arc  increased  the 
variability  of  M2  conversion  at  the  Luzon  Strait.  We 
estimate  the  westward  energy  flux  from  the  Mariana 
Arc  as  1.3  kW  m-1  from  Kerry  et  al.  (2013),  assuming  a 
line  source  and  no  dissipation  or  spreading.  This  is 
—4%  of  the  mean  semidiurnal  coherent  flux  at  Al  and 
is  likely  an  upper  bound.  In  the  absence  of  detailed 
stratification  or  velocity  information  along  the  entire 
transect,  we  assume  that  the  phase  speed  and  travel 
time  vary  similarly  to  a  LZS  model  point  east  of  the 
ridges.  The  phase  of  the  arriving  wave  has  a  range  of 
more  than  360°,  and  the  coherence  of  the  semidiurnal 
flux  in  this  2D  model  is  69%.  This  is  lower  than  the 
observed  semidiurnal  coherence  at  both  Al  (93%) 
and  S9  (73%),  so  this  mechanism  could  account  for 
those  values.  The  coherence  of  semidiurnal  conver¬ 
sion  at  Al  in  this  model  is  86%,  close  to  the  observed 
value  of  88%. 

Though  these  simple  2D  calculations  neglect  possible 
feedbacks  that  could  change  conversion,  and  do  not  in¬ 
clude  any  scattering  or  dissipation,  they  provide  a  rough 


estimate  of  the  possible  magnitude  of  incoherence  that 
these  mechanisms  could  produce.  The  results  suggest 
that  the  arrival  of  remote  semidiurnal  internal  tides 
from  the  Mariana  Arc  could  account  for  the  observed 
semidiurnal  incoherence  at  both  Al  and  S9  and  that  the 
reflection  of  the  diurnal  internal  tide  from  the  conti¬ 
nental  shelf  could  account  for  the  observed  level  of  di¬ 
urnal  incoherence  at  S9.  However,  the  observed  —30% 
coherence  of  the  diurnal  internal  tide  at  Al  cannot  be 
explained  by  these  simple  2D  models. 

To  investigate  further,  we  use  the  3D  LZS  model  and 
compute  coherence  in  the  same  way  as  for  the  obser¬ 
vations  (Fig.  16).  Coherence  in  both  frequencies  is  lower 
along  the  eastern  ridge,  especially  south  of  20°N.  The 
diurnal  band  (Fig.  16c)  also  has  a  region  of  low  co¬ 
herence  extending  to  the  northwest  from  the  eastern 
ridge  near  20.5°N.  Both  diurnal  and  semidiurnal  internal 
tides  are  largely  coherent  in  the  area  of  S9.  Al  is  located 
in  a  region  of  low  diurnal  coherence  and  higher  semi¬ 
diurnal  coherence.  A  zonal  transect  (Fig.  16a)  through 
Al  at  20.5°N  shows  excellent  agreement  between  model 
and  observations.  A  transect  through  S9  (Fig.  16b) 
shows  good  agreement  in  the  diurnal  band  and  slightly 
higher  than  observed  coherence  in  the  semidiurnal 
band.  Since  the  model  reproduces  the  observed  pat¬ 
terns,  we  can  use  it  to  explore  mechanisms  for  lowering 
coherence.  Notably,  since  the  model  does  not  include 
remote  waves  from  the  Mariana  Arc,  more  local  effects 
appear  responsible. 

A  clue  to  the  diurnal  incoherence  at  Al  is  given  in 
Buijsman  et  al.  (2014),  who  identified  a  large  gradient  in 
the  phase  of  mode-1  diurnal  pressure  at  that  location. 
We  hypothesize  that  incoherence  arises  as  the  variable 
Kuroshio  and  mesoscale  shifts  this  pattern  back  and 
forth,  which  Buijsman  et  al.  (2014)  were  not  able  to 
examine  since  the  MITgcm  model  did  not  contain  me¬ 
soscale  forcing.  Here,  we  compute  the  phase  of  mode-1 
perturbation  pressure  from  the  LZS  model  using  a 
sliding  harmonic  fit  to  the  diurnal  and  semidiurnal  fre¬ 
quencies  in  overlapping  3-day  windows.  The  diurnal 
phase  (Figs.  17b, e)  has  a  large  gradient  that  runs  from 
northwest  to  southeast,  passing  near  Al  before  con¬ 
tinuing  south  along  the  eastern  ridge.  This  pattern  and 
strong  north-south  gradient  between  the  ridges  is  simi¬ 
lar  to  that  seen  in  Buijsman  et  al.  (2014)  but  appears  to 
be  less  zonal  near  20.5°N,  possibly  because  of  the  addi¬ 
tion  of  the  Kuroshio  in  the  LZS  model.  The  gradient  of 
the  semidiurnal  phase  (Figs.  17c, f)  is  smaller  at  Al,  in 
agreement  with  the  higher  coherence  of  that  compo¬ 
nent.  Near  the  location  of  S9,  the  gradient  of  diurnal 
phase  is  small  and  the  semidiurnal  phase  gradient  is 
larger,  consistent  with  the  relative  levels  of  coherence 
observed  there. 
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Fig.  16.  Percent  variance  of  the  mode-1  energy  flux  explained  by  coherent  component  in  LZS 
model  for  the  (c)  diurnal  and  (d)  semidiurnal  frequency  bands,  (top)  Cross  sections  through 
(a)  A1  and  (b)  S9,  (bottom)  with  location  indicated  by  dashed  lines.  Observed  values  from 
moorings  are  plotted  as  squares.  Bathymetry  along  cross  sections  is  shown  in  gray,  with  a  scale 
at  lower  right  of  (b). 


To  test  our  hypothesis  that  shifting  of  this  gradient  is 
responsible  for  the  lower  coherence,  we  also  examined 
the  standard  deviation  of  the  baroclinic  phase  (not 
shown).  It  was  generally  largest  where  the  largest  gra¬ 
dients  exist,  consistent  with  this  spatial  pattern  being 
advected  or  shifted  laterally. 

The  shifting  of  the  internal  tide  phase  pattern  must 
be  caused  by  changes  in  stratification  or  velocity  shear, 
which  in  the  Luzon  Strait  are  dominated  by  the  in¬ 
trusion  of  the  Kuroshio.  Figure  17  shows  the  low- 
frequency  velocity  at  100  m  at  two  different  times 
(Figs.  17a, d).  In  the  first  period  (top  row),  there  is  a 
weaker,  more  diffuse  Kuroshio,  and  the  white  diurnal 
phase  contour  passes  through  A1  (Fig.  17b).  In  the 
second  time  period  (bottom  row),  the  Kuroshio  is 
stronger  and  more  jetlike,  and  the  diurnal  phase  pat¬ 
tern  around  A1  (Fig.  17e)  appears  to  be  shifted  to  the 
north.  This  suggests  that  the  shifting  of  the  diurnal 
pressure  phase  pattern  is  related  to  the  low-frequency 
velocity  or  Kuroshio.  However,  the  specifics  of  the 
interaction  between  the  3D  internal  tide  response  and 


the  background  flows  are  complicated  and  beyond  the 
scope  of  this  paper. 

6.  Conclusions 

The  coherence  and  variability  of  internal  tides  in  the 
Luzon  Strait  were  studied  using  measurements  from  two 
moorings  deployed  for  50  days  during  the  IWISE  ex¬ 
periment.  The  observations  reveal  an  extremely  ener¬ 
getic  internal  wave  field,  with  depth-integrated  diurnal 
(semidiurnal)  energy  fluxes  up  to  54  (63)kWm-1.  Ve¬ 
locity  spectra  measured  at  the  moorings  were  elevated 
up  to  15  times  above  the  theoretical  GM  spectrum 
(Fig.  3),  likely  because  of  energy  input  from  strong  in¬ 
ternal  tides  and  implying  elevated  turbulent  mixing 
rates.  The  modal  distribution  of  internal  tide  energy  flux 
was  red,  with  most  of  the  flux  contained  in  the  first 
two  modes. 

The  temporal  variability  of  the  internal  tide  is 
quantified  by  computing  the  coherent  component, 
described  by  a  sum  of  tidal  frequencies  with  constant 
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Fig.  17.  (left)  Low-frequency  velocity  at  100  m  and  phase  of  (middle)  diurnal  and  (right)  semidiurnal  mode-1 
pressure  for  (top  and  bottom)  two  time  periods  of  the  LZS  model.  The  same  phase  value  is  contoured  in  white  to 
show  the  shifting  of  the  phase  pattern.  A  constant  offset  is  added  to  phases  in  the  second  time  period  to  account  for 
different  barotropic  phases. 


amplitude  and  phase  over  the  50-day  mooring  de¬ 
ployment.  The  semidiurnal  internal  tide  at  A1  and 
both  frequencies  of  internal  tide  at  S9  were  largely 
coherent  (Table  2),  with  more  than  70%  percent  of  the 
energy  flux  variance  explained  by  the  coherent  com¬ 
ponent.  In  contrast,  the  diurnal  internal  tide  at  A1  is 
significantly  less  coherent,  with  the  coherent  energy 
flux  only  30%  of  the  total.  The  phase  and  amplitude  of 
diurnal  perturbation  pressure  at  A1  varies  significantly 
over  the  mooring  deployment,  resulting  in  the  low 
coherence. 

Simple  2D  calculations  suggest  that  reflection  of  the 
diurnal  internal  tide  from  the  continental  shelf  and 
remote  semidiurnal  internal  tides  arriving  from  the 
Mariana  Arc  could  account  for  the  observed  levels  of 
coherence,  with  the  exception  of  the  diurnal  internal 
tide  at  Al.  Instead,  the  much  lower  coherence  of  the 
diurnal  internal  tide  at  Al  is  a  result  of  the  mooring’s 
location  near  a  large  gradient  in  the  3D  pattern  of  the 
baroclinic  pressure  phase.  The  LZS  model,  which 
contains  a  realistic  Kuroshio  and  mesoscale  variabil¬ 
ity,  shows  that  changes  in  the  Kuroshio  intrusion  into 
the  Luzon  Strait  occur  at  the  same  time  as  shifts  of  the 
internal  tide  phase  pattern  in  the  Luzon  Strait. 


However,  the  multiple  sources  and  3D  nature  of  the 
interference  pattern  makes  it  difficult  to  isolate  the 
details  of  this  mechanism.  Continued  progress  will 
require  a  combination  of  data  from  many  different 
observation  platforms,  combined  with  models  such  as 
the  LZS  that  include  the  Kuroshio  and  mesoscale 
variability. 

Acknowledgments.  IWISE  was  funded  by  the  Office 
of  Naval  Research  (ONR).  We  thank  the  Captain  and 
crew  of  the  R/V  Revelle  and  R/V  OR  1  for  their  skill  and 
hard  work  to  acquire  these  measurements.  We  would 
also  like  to  thank  John  Mickett,  Eric  Boget,  Zoe 
Parsons,  Paul  Aguilar,  Tom  Peacock,  Amy  Waterhouse, 
Ruth  Musgrave,  Hayley  Dosser,  Ke-shien  Fu,  and 
Chungwei  Lu  for  their  work  in  deploying  the  moorings 
and  making  shipboard  measurements.  This  work 
was  supported  by  ONR  Grants  N00014-09-1-0219 
(Dr.  Alford  and  Rainville),  N00014-1 1-1-0184  (Pickering), 
N00014-05WX-2-0647  (Dr.  Ko),  ONRDC32025354 
(Dr.  Buijsman),  and  N00014-09-1-0281  (Dr.  Nash).  Com¬ 
ments  from  two  anonymous  reviewers  helped  improve 
the  manuscript.  IWISE  would  not  have  been  possible 
without  the  hard  work,  talent,  and  generosity  of  our 


June  2015 


PICKERING  ET  AL. 


1593 


Fig.  Al.  Figure  illustrating  the  removal  of  mooring  blowdown 
from  E-chain  measurements  at  S9  and  the  effect  on  calculated  ver¬ 
tical  displacement,  (a)  Zonal  barotropic  velocity,  (b)  Thermistor 
string  temperature  plotted  vs  nominal  depth  of  instruments.  Black 
lines  indicate  depth  of  SBE-37  and  ADCP  at  top  and  bottom  of  the  T 
chain,  respectively,  (c)  Temperature  corrected  for  mooring  blow¬ 
down  using  pressure  measured  at  top  and  bottom  of  the  T  chain.  The 
15°C  isotherm  is  contoured  in  red  in  (b)  and  (c)  to  illustrate  the  effect 
of  blowdowns  on  inferred  vertical  displacement. 
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APPENDIX 

Estimating  Internal  Wave  Displacement  from 
Temperature  Chain  Measurements  with  Mooring 
Knockdown 

Extremely  strong  upper-ocean  currents  in  the 
Luzon  Strait  inhibited  moored  profilers  from  climb¬ 
ing  at  times  during  the  pilot  experiment.  To  prevent 
this  issue  during  the  main  experiment,  mooring  S9 
was  redesigned  with  an  upward-looking  ADCP  and  a 
temperature  chain  to  measure  the  upper  300  m.  The  T 
chain  consisted  of  SBE-56  temperature  loggers 
spaced  about  15  m  apart,  with  an  upward-looking 
ADCP  at  the  bottom  of  the  T  chain  and  a  SBE-37 
MicroCAT  at  the  top  of  the  T  chain.  Here,  we  discuss 
the  correction  of  mooring  knockdown  in  these 
measurements. 


Mooring  knockdown,  determined  from  pressure 
measured  at  the  top  (SBE-37)  and  bottom  (ADCP)  of 
the  T  chain,  was  up  to  —75  m  during  strong  flows  at 
spring  tide.  The  instantaneous  depth  of  each  SBE-56 
temperature  sensor  was  calculated  from  the  pressure 
recorded  at  the  top  and  bottom  of  the  T  chain,  assuming 
that  the  T  chain  remained  a  straight  line  between  the  top 
and  bottom.  This  assumption  can  be  tested  by  compar¬ 
ing  the  difference  in  depth  between  the  top  and  bottom 
(5z)  to  the  known  length  of  the  wire  L  =  237  m.  The 
maximum  value  of  (L  -  Sz)  was  about  10  m  or  about  4% 
of  L,  indicating  that  this  is  a  good  assumption.  The  cable 
could  be  slightly  curved  instead  of  in  a  straight  line. 

Figure  Al  shows  temperature  plotted  at  the  nominal 
depth  of  the  T-chain  instruments  (Fig.  Ala)  and  after 
interpolating  to  the  actual  depth  (Fig.  Alb).  Without 
correction,  large  mooring  blowdowns  during  spring  tides 
(i.e.,  day  183)  would  be  interpreted  as  extremely  large 
vertical  displacements.  These  vertical  displacements 
have  a  large  effect  on  the  internal  tide  energy  flux. 
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